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the upper critical dimension, they can when dimensions are lowered. 

PACS numbers: 64.60.Ak, 64.60.Fr, 75.10.Nr 



I 

o 
o 



> 

O 
m 
o 



X3 
O 

o 



X 



I. INTRODUCTION 

The random-field Ising model (RFIM) is the Ising 
model coupled to a random magnetic fieldQ. Although 
it has been extensively studied for about three decades, 
there remain many unresolved problems 0, one of which 
concerns dimensional reduction!^. 

Dimensional reduction claims that the critical behav- 
ior of the d dimensional RFIM is equivalent to the d — 2 
dimensional pure Ising model. Following this argument, 
phase transition does not occur in the three-dimensional 
RFIM. However, it is rigorously proved that the phase 
transition does occur in this model, so that dimensional 
reduction does not hold at least d = 30. Various numer- 
ical computations are performed to obtain critical expo- 
nents in three dimensions [^0,0, II]- On the other hand, 
it is not understood whether dimensional reduction works 
in other dimensions lower than the upper critical dimen- 
sion, which is believed to be six. 

Since standard perturbation for the Ginzburg-Landau 
(GL) theory of the RFIM simply leads to the result 
consistent with dimensional reduction, other approaches 
were applied. Schwartz et al proposed modification of di- 
mensional reduction, which indicates correspondence be- 
tween the d dimensional RFIM and the pure Ising model 
in d' = d—2 — ri{d') . Mezard and Young suggested the 
replica symmetry breaking of the RFIM by extrapolation 
from 1/m expansion, where m is the components of the 
spin lOj. Lancaster et al computed exponents, pa ying 
attention to many solutions of mean field equations [ij- 
Although these works support the breakdown of dimen- 
sional reduction in other dimensions, it has not yet been 
settled. 

The breakdown of standard perturbation may be 



caused by overlooked relevant operators. In 4 + e dimen- 
sions. Fisher |l3) and Feldmanfl^ pointed out that there 
are infinitely many relevant operators in n-component 
random spin models. Near the upper critical dimension, 
Brezin and de Dominicis investigated that the GL Hamil- 
tonian corresponding to the RFIM has not only the usual 
(j)'^ interaction 1231 
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where a denotes the replica index, but also the following 
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It means that the space of the coupling constants extends 
to the five dimensions and that a renormalization-group 
(RG) trajectory should be considered in five-dimensional 
space. They claimed that dimensional reduction does 
not work even near the upper critical dimension since 
the non-trivial fixed point of 0(e) becomes unstable in 
d = 6-e. 

However, their analysis has ambiguity that originates 
from the zero-replica limit. In fact, they proposed two 
different limiting procedures and derived two sets of beta 
functions quantitatively different each other. 

In this paper, we circumvent the ambiguity and recon- 
sider RG flow in the extended coupling-constant space. 
To this end, we use RG for the random probability dis- 
tribution that controls random potentials including the 
magnetic field. This method was initiated by Harris and 
Lubensky 1^, il8„ .19] in order to investigate a randomly 
diluted magnetic system with non-magnetic impurities. 
In this case, it is confirmed that the replica method 
is consistent with the Harris-Lubensky method J,8, 2Cj. 
Here we extend the random probability distribution 
adopted in the previous literatures^, 0, \m. It corre- 
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spends to taking into account the replica interactions (O 
as well as 

According to Refs.0,0, the couphng constants in 
the Hamiltonian are extended to inhomogeneous ran- 
dom potentials that are correlated to each other following 
some random probability distribution P. After perform- 
ing standard RG in the inhomogeneous potentials, the 
change of the Hamiltonian can be pushed into transfor- 
mation of P. Distribution P is characterized by non- 
trivial cumulants with some parameters. Therefore the 
change of P defines RG flow in that parameter space. We 
can investigate critical behavior from RG flow. An ad- 
vantage of this method compared to the replica method 
is that any ambiguities originating from limiting proce- 
dures do not arise. 

This paper is organized as follows: we show in the 
next section the RG for probability distribution of the 
Gaussian random magnetic impurities. We can naturally 
introduce five parameters that completely correspond to 
the five coupling constants for the interactions 1^ and 
^ in the replica Hamiltonian. In Sec. IIHI we ana- 
lyze the RG by perturbation and derive the recursion 
equations of the parameters that specify random-field 
distribution. We also compare the result with that of 
the replica methodjl^]. In Sec. IIVI we introduce ex- 
pansion parameters of physical quantities that are con- 
venient for analysis near the upper critical dimension. In 
the language of the replica method, it means to deter- 
mine the scaling dimensions of the operators in and 
(O. We show that the interaction is the unique rel- 
evant operator near the upper critical dimension in all 
the 0^ interactions. Then arguments of dimensional re- 
duction can survive near the upper critical dimension. 
However, as the dimensions are lowered, other interac- 
tions can become relevant, which cause the breakdown 
of the dimensional reduction. This picture is consistent 
with hig h-temperature expansion studied by Houghton 
et al^^. The last section is devoted to summary and 
discussion. 



II. RENORMALIZATION-GROUP 
TRANSFORMATION FOR THE RANDOM 
PROBABILITY DISTRIBUTION 



In this section, we formulate the renormalization-group 
transformation (RGT) for the random probability distri- 
bution in the presence of magnetic impurities following 
Refs.|Tl[l3. The GL Hamiltonian to the RFIM is usu- 
ally described by 0: 



l3H^l-f {kl+t)(l>(h)cf>{-h)+ f vi{ki)^{ki) 

^ Jki Jki 



Ui 

'4! 



<^(fcl)</'(fc2)0(fc3)0(-/Cl-fe-fe) 



fel ,fc2 ,fc3 



where the momenta fci , ^2 , ^3 belong to 
IC = {k : < |fc| < A} 
and the integral means 

d'^ki r (fiki 



f = f ^zl... f 

Jki, ..Mi Jkieic (271')'' JkieK. 



(4) 



(5) 



The random magnetic field vi (fc) follows the Gaussian 
distribution Po[t'i] 



Vvi = Y[ dvi{k). 



(6) 



fce/c 



We treat the random field vi as an external field for a 
while and examine the RGT. Namely, we first integrate 
higher-momentum components of cf) and then perform the 
rescaling of the potentials and the fields appropriately. 
Let us observe correction terms that appear in higher- 
momentum integration. Let G{q) be the free propagator: 
G{q) = {q^ -\-t)~^. By integrating the higher- momentum 
components by perturbation in ui, we have the following 
correction to the 0^ term in the leading order; 



Ui j Vi{q)G{q)vi{pi+p2 + q)G{pi+p2 + q)<t>{pi)(t>{p2)- 

Jq 

(7) 

Here integration is carried out on higher-momentum 
space /C> and 4){pi) is a lower-momentum component. 
Similarly, the 4>'^ term gets 
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G{pl+P2+ql)Y[<l>{p^)■ (8) 



k=l 



Further, a cubic interaction emerges through, for exam- 
ple, 

> / 3 3 \ 3 



G{pi +P2 + qi 



(9) 



k=l 



(3) 



These corrections are depicted in FIG. Q Eqs. Q), © 
and @ respectively correspond to (a), (b) and (c) of 
the figure. The correction terms indicate that the coef- 
ficients of the (jP and 0* terms in the Hamiltonian no 
longer possess translational invariance. In order to ab- 
sorb those terms, we need to extend (fc^ -I- t) and ui to 
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(a) (b) 




(c) 



FIG. 1: Graphical representations for Eqs. Q, ||HJ and 0. 
A solid line stands for G{q) and a dashed line for <j){Pi)- The 
higher-momentum components of the random field vi is de- 
picted by an open circle. 

inhomogeneous potentials. In addition, the cubic interac- 
tion should be also included. Thus the GL Hamiltonian 
S{= f3H) is generalized as 

5['^;p]=En/ viih,...,ki)(bih)---cl,{ki), (10) 

1=1 Jki,...,ki 

where p represents all of the inhomogeneous potentials: 

P = {vi,V2,Vz,V4) . (11) 

Starting with the Hamiltonian (|10|l . we describe the 
RGT definitely. The higher and lower-momentum spaces 
are respectively defined as 

K.> = {q:L-^K<\q\<K] (12) 

and 

/C< - {p : < IpI < L^^A} (13) 

with L > 1. The Hamiltonian S is decomposed into 
and 5^, where 5^ consists only of the lower momentum 
components of the field. The remaining part is denoted 
by 5>,i.e., 

S = S<+S>. (14) 

Let us integrate over the higher-momentum components 
(/)(g) in the partition function. Namely, 




where 5S can be written as 



6S = -\n Yl d<j){k)e~'^^ 
1=1 ■ '' 

-\- irrelevant terms. (16) 
Next the scaling transformation is performed as 



</.(p)=LV'(fc), (17) 



where k is related to p by 



k = Lp. (18) 



Defining the new inhomogeneous potentials v'^ by 



vl (ki, ki) = L'C-'^) {in (pi, ...,pi) + Svi {pi, ...,pi)) , 

(19) 

the Hamiltonian turns back to the original form: 



s< + ss ^ [ "'(^1' fcO-^'lfci) • ■ • 

1 = 1 Jki,...,ki 

= S[c^';p'], (20) 



where p' — (wj^, ^3, W4). We then get the RGT for the 
potential p ^ p' . Furthermore, a correlation function 
transforms as 



(0(pi), • ■ • , <j){Pn)) S[^:p] = ^ (-^(fcl), ■ ■ ■ , ^{kn)) S[d,;p'] : 

(21) 

where means the thermal average using the Boltz- 
mann weight . 

We proceed to the average over the random potentials. 
It should be noted that Eqs. iQ, © and © show that 
V2, V3 and W4 come to have non-trivial correlations with 
each other. Therefore we can regard them as random po- 
tentials that obey some probability distribution as well as 
fi. We shall denote probability distribution by P[p]. For 
example, if we take Eq. as the initial Hamiltonian, 
the corresponding distribution Pi is given by 
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Pl[p] = Po[vi] H 5[v2ih,k2)- {kl+t){27rfSih+k2))Siv3ih,k2,k3)) 



fci ,k2 ,fc3,fc4 



6 V4 (fci, fc2, fe, ^4) - ui {2n f 6 kj 



(22) 



The RGT for the random potential p ^ p' can be 
pushed into change of P by the foUowing rule: 



added to Eq. 



J F[p']P[p]Vp = J F[p]P'[p]Vp, 



(23) 



where F is an arbitrary functional of p and 



^?p=n n dviiki,...,ki). 



Namely, P' is formally written asp^ 

dp{p') 



(24) 



py] = p[p{p') 



dp' 



(25) 



Eq. (|23J) defines RGT P ^ P' , keeping the Hamiltonian 
invariant. For instance, using Eq. (|21|) . the random- 
potential average of the correlation function becomes 



P[P] {4>{Pl), • ■ ■ , <l>{Pn)) S[4,;p] 

P'[p]L"'{Hk^),.-- ,<l){kn))s[^.^p]'Dp. (26) 



In the practical computation of the random-potential 
average, we give cumulants among random potentials in- 
stead of giving explicit form to the probability distribu- 
tion. First, as in the original theory, we put 

[vi{k,)-v,{k2)]p^p-^ = A {27rf S{k, + k2) 
[f2(A:i, /C2)]p[^] = {kf + t) {2t:Y d{ki + k2) 

[z;4(A:i,fc2,fc3,fc4)]p[^] =ui {2T:f 5 I j '(2"^) 



where [ ■ Jp^^j means to take the average over the random 
variables with distribution P[p]. The semicolon in the 
bracket means the cumulant product: e.g., [X;y]pj^j = 
[XFjpj^j — [-'^Ipjp] [^]p[p]- We assume that there is no 
long-range correlation between random variables. In this 
assumption, the following non-vanishing cumulants are 



/ bl(^l);W3(fc2,fc3,^4)]p[p] \ 

'V2{ki,k2);v2{k:i,k4)]p^p^ 
'v2{ki,k2); viiks); vi(fc4)]p[p] 
\ [vi{ki);vi{k2);vi{k3);vi{k4)]p[p] J 




{27rfd 



Note that we respect the symmetry 




Vl 



-wi, 0- 



so that 



Mkl)]p[p] = b3(^l,^2,fc3)]pf^, =0 



Plp] 



(28) 



(29) 



(30) 



The parameters Uj{j — 1, ...,5) have the scaling dimen- 
sion A — d (measured by the inverse length) . Cumulants 
other than Eqs. H27|) and H28|) are ignored since they 
are associated with higher-dimensional parameters[2^. 
Probability distribution P is characterized by parameters 
(t. A, ui, U5) = p. Transformed distribution P' is sim- 
ilarly characterized by (i'. A', u'j^, Ug) = p' . Namely, p' 
is obtained taking the random-potential average of Eqs. 
(|77jl and (EHl) in use of P'[p\ instead of P[p\. 

Transformation p 1-^ p' is analyzed using the following 
formula: 



hi(^i\ 



' ^/i); ' ' ■ ; '^in (^"' )] p[p] ' ('^-'-) 



which is easily checked from Eq. H23|) . We apply it to the 
non-vanishing cumulants. By definition, the left-hand 
side can be written in p' . On the other hand, if the 
transformed potentials are described by the original ones 
employing Eq. H19|l . the right-hand side can be evaluated 
in terms of p. In this way, p' is expressed in terms of p. 



III. PERTURBATIVE RENORMALIZATION 
GROUP 

A. Diagrammatic expansion 

In this section we obtain transformation ^ 1— > by 
perturbation. We consider all UjS as small parameters 
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and express ^' up to 0{uiUj). On the other hand, we 
do not regard t and A as small since they are relevant 
parameters of the dimension 2. 

We can formulate the perturbative expansion by 
adopting the random-potential average of 112(91,92) and 
the linear term as the unperturbed Hamiltonian 5*0: 

So = I [ [«2(fcl,fc2)]p[p]0(fcl)<^(fc2)+ / V,{-k)m 
^ Jki,k2 Jk 



^ \k^ + t)^{k)^{-k)-\ / vi{~k)G{k)vi{k), 



(32) 



where 



(^(fc) = (^{k) + G{k)vi{k) 



(33) 



The remaining terms are denoted by V and treated as 
the perturbative Hamiltonian. Let 



J=2 



(34) 



where 



^ JkiM ^ ' 

Vj = ^ Vj{ki,...,kj)(j){ki)---(l){kj), 

J ' J ki ,. . . .kj 

(j=3,4). (35) 

As we have performed in Eq. 114|l . we divide Sq and 
V into 



So — Sq + Sq 

V = V> + V< 



(36) 



respectively. The free part of the partition function is 
denoted by 



Zo = J Y[ dc^{q) 

qeK> 



^0 



Writing 



(37) 



(38) 



6S in Eq. 1)16(1 is expanded as 

-SS ^ logZo + log(e'^' 



\ogZQ-{V>)^ + -{V>;V>)^^ 



3! 



{V>;V>;V>)^+0{V>\ (39) 



V2 + 



FIG. 2: Corrected potential 112. An internal line carrying 
higher momentum is represented by a solid line, while an ex- 
ternal one by a dashed line. An open circle stands for vi{q). 






FIG. 3: four-point vertices 



where the semicolon represents the cumulant with respect 
to the thermal average ( )^. Note that, in practical com- 
putation, (j){q) in is understood as ip{q) — G{q)vi{q) 
and the Wick theorem is applied to ip{q). Comparing 
Eas. H16|l and ((39(1 . we can obtain the perturbative ex- 
pansion of 6vi{pi, ...,pi), which leads to transformed po- 
tential v'l defined in Eq. ((19(1 . Inserting the resultant v'^ 
into the right-hand side of Eq. (|^ , we can express the 
transformed parameters in terms of the original ones. In 
order to proceed with this program efhciently, we will use 
the diagrammatic technique. 

Let us depict v'l by a shaded circle with I dashed lines. 
FIG. 121 exemplifies the case of ^ = 2. The shaded circle 
contains Svi, which is expanded by connected diagrams, 
according to Eq. 1(39(1 . The random potential i;;(fci, fc;) 
makes an Z-point vertex. In particular, vi{q){q € /C>) is 
expressed by an open circle. It should be noted that total 
momentum is not conserved at the vertex vi{ki, ki). 

When we take the random-potential average, non- 
vanishing cumulants [w4]p[p], b3;wi]p[p], b2;w2]p[p], 
[w2;wi;wi]p[p] and [vi;vi;vi;vi]p^p^ can be graphically 
represented as four-point vertices at which the total mo- 
menta, carried by the four lines, are conserved due to the 
delta function in Eq. See FIG. El 

Further, the cumulant [wi;wi]p[p] tells us that two in- 
ternal lines ended with the open circles are merged by 
the random-potential average and produce A. We here 
depict it by a cross, x , as in FIG.^ The cross reminds us 
that the two internal lines are not connected with respect 
to the thermal average. 
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many internal lines as possible. After we obtain the set 
of all one-loop diagrams, we have to compute their nu- 
merical factors. Although the procedure is fulfilled in 
appendix here we present the following examples. 



FIG. 4: Merging two open circles 



B. Example 



The right-hand side of Eq. is graphically repre- 

sented as the left-hand side of FIG. |5l When we expand 
the shaded circles in FIG. El we get a summation of the 
diagrams that contain several connected components, as 
shown in the right-hand side of the figure. If we want 



Let us consider FIG. Ela). Putting x on the inter- 
nal Hues, (b), (c), and (d) in the figure are generated. 
However, (c) and (d) are excluded because they contain 
connected components without external lines. Counting 
the external lines of the connected components in FIG. 
|SJa), we find that it makes a contribution to u'^, while 
FIG. El^b) contributes to Ug. First we compute the con- 



= V + :<:>:' + ^:<y:' + 



(a) 



(b) 



0*. 



= yt + ;<:x + + 



FIG. 5: Perturbative expansion for cumulants contributing 
^ (i = 1, ■••,5) 

to know in 0(ujUk), we need to consider the one-loop 
diagrams made of the vertices in FIG. 13 and x, includ- 
ing the numerical factors of those diagrams. Note that a 
one-loop diagram that contains a connected component 
without external legs does not appear in the right-hand 
side of FIG. El because the shaded circles in the left-hand 
side consist of only connected diagrams. In particular, it 
is ruled out that there are two or more x s on an internal 
line. 

In order to write down all one-loop diagrams in the 
right-hand side of FIG. |S1 we must first consider all the 
admissible diagrams without x , and then put x on as 



(c) (d) 
FIG. 6: Diagrams generated from (a) by putting x. 

tribution from FIG.|ni[a) to 

u', i2nf 6 (^^k}j = [v',ik,,k,),v[ik,y,v[ih)]p^^yiAO) 
Using Eq. ifT^ll . FIG.E^a) comes from 

(41) 

The second term arises because the roles of v'l^k^) and 
v'i{k4) can interchange. The two terms in the right-hand 
side of Eq. H41I) equally contribute to U4. In general, if a 
cumulant contains multiple potentials of the same kind, 
the number of ways of interchanging should be counted. 
We shall denote the number with nn. As for the present 
case, 



nji = 2. 

Now we compute the first term. Defining 



(42) 



V, 



(m) 



1 / I 



{pi,...,Pm) = - I I X 



l — ra 

/ VI {pi, ...,pm,qi, ..,qi-m) Yl 7(43) 

Jqi,...,qi-^ 
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we show explicitly the terms in Sv2{pi,P2) necessary for 
making FIG.IHIa) as 

Sv2ipup2) - -2\-^{v^'\piyM'\p2))^ + --- 



ai / V2ipi,qi)v2ip2,~qi)G{qi) 

Jqi 



(44) 



Here the dots in the right-hand side indicate terms not 
employed for FIG. Ela). In the first equality, the factor 
2! originates from the normalization of Sv2{pi,P2), while 
1/2! comes from the expansion of e~^^. The minus sign 
is so put because it is a second-order perturbation term. 
Performing the thermal average, we can easily check that 



ai 



(45) 



Similarly, 

Svi{p) 



0-2 / V2{p3,q2)vi{-q2)G{q2) 
•Jq2 



(46) 



with 



0-2 



(47) 



The factors ai and 02 are associated with the connected 
components of the diagram we have computed. We shall 
denote the product of them with ua, i.e.. 



Ua = 0x02 = 1. 



(48) 



Next we go to the random-potential average. Inserting 
Eqs.lI33|) and (gEl into the first term of Eq. the 
random-potential average becomes 



-i{e-d) 



G{qi)G{q2) [^2(^1, gi)w2(p2, -Qi); ^2(^3, 92)wi(-<?2); "1(^4)] 



91,92 
> 



P{p\ 
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G{qi)G{q2) [^2(^1, 9i); ^2(^3, 92)] p[p] [^2(^2, -gi); wi(-g2); ^'1(^4)] p[p] 



+ {v2 {pi,qi) ^ V2 {p2, -qi)) 



Tie-3d 

= L 



[G{q)G{p2 +PA-q) + G{q)G[pi +Pi + q)) 
I 



(49) 



Note that there are two ways of contraction with respect 
to the random-potential average in the first equality. We 
denote the number of ways of the contraction by nc- 
That is, 



nc = 2 



(50) 



in the case of FIG. Efa). 

Thus the correction to U4 by way of FIG. Ela) is read 
from Eq. (|49|l . combining with the second term of Eq. 
(glj. Letting p.j = (z = 1, ...,4) in G, we find that the 
numerical factor n p associated with the diagram is 



np = nATiRnc = 4. 
Therefore the correction to U4 is 



(51) 



4e-3d 



U3U4 I G {q). 

(52) 



Let us turn to FIG. Elb), which affects cx 
[v[{pi);v[{p2);v'i{p3);v[{p4)]p^py For creating this dia- 
gram, we need three ^V'2*'^^(p)^ and one vi(j>). Obvi- 
ously, 



riR = 4. 



(53) 



The numerical factor associated with the connected com- 
ponent was already computed in Eq. H4fi|l . Their product 
becomes 



nA = (-1) X 1 



1. 



(54) 



Since the number of contractions is equal to the number 
of ways of contracting V2 (p, q) and vi (p) , 



nc 



(55) 
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Thus the contribution to Ug of FIG.|Hfb) is 



npAu2U3L 



(56) 



equation in 0{uiUj). The coefficients ua, nc, and nR 
associated with each diagrams are presented in appendix 
IXI Here we present the result. 



with 



np = nRncriA 



Results 



-12. 



(57) 



u'^ = L^'-'" {u, + Su,) , (j = l,...,5) 
t' = L^o-'i {t + St) 
A' = L^e-d (A + 6 A) 



(58) 



After writing all one-loop diagrams appearing in the 
right-hand side of FIG. [SJ we can obtain the recursion 



with 



Sui 

SU2 
SU5 



3A3A + ^A2 ] u\ -I- 6^2W1W3 



-'iAzAu\U'i — -A2U1U2 + 6A2U2U2 — 3A2U2 + SA2U1U4 



-^4A u-^ — (2^3A + A2) uiu^ + A2U1U2 + A2U1U4 + 4^2^*3 ~ 3A2U2 — 6A2U2U3 

2 f 1 \ 2 1 

AaA U1U3 - ( S^sA -I- -A2 I U1U4 - 4^3A{t3 + -A2U1U5 - -42M2U3 + 8^2lt3"4 

6^4A^U3 — 36^3AU3U4 — 6A2U2U4 + 6^2^31*5 -I- 24:A2u'l, 



(59) 



where we have defined 

An 

and 



GiqT 



(60) 



Similarly, 5t and 5 A are given as 



U3=U2+ U3. 



(61) 



St = ^{A2A + Ai)ui~AiU3 



SA 



:^2^3 



:B2 A' 



tAiAs 



Bi + -Ai]A- -A1A2 - -60 ) < 



+^AiA2UiU2 + {{A1A3 + Al + 2Bi) A + A1A2 + Bo) mus 

- {A1A2 + Bo) U1U4 - {A1A2 + Bo) ul 

A2AU3 + U2A1 - 2U4A1 + ^B^A^ul - {]^A^ + \AiA2 + ^So^ U1U2 

- {{A2A3 + 2S2) + (^1^3 + Bi) A) U1U3 

+ {{Al + 3Bi) A + A1A2 + Bo) U1U4 - \bouiU5 - Bouj 

+ {2A1A2 + 2Bo) U2U3 + (3i3i + ^2 + 2^1 A3) Au^ - (4^i ^ + 4So) W3M4, 



(62) 



where 



Bo = 
Bi = 



Giqi)G{q2)Giqi + q2) 



91,92 

> 



D. Comparison to the replica method 

We compare the recursion relations shown above to 
those obtained by the rephca method0,|22, where the 



G(qi)^G{q2)G{qi+q2) 
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perturbative Hamiltonian is 

rep ^1 rep 

rep rep 

+\- E ^l^l + ^ E 

a,/3— 1 a,/3,7— 1 

rep ''^ 

+ ^ E <^"<^/3'/'7'^5- (64) 
Q,/3,7,(5— 1 

Brezin and de Dominicis proposed the two limiting pro- 
cedures listed below: 

(A) u'l"^, Ug"'' are fixed, then take n ^ 0. 

(B) u^^'^ ,nu^2^ ,nu^^^ ,n^u^^ ,n^u^^^ are fixed and take 

We can check that the recursion equations H58() - H61|l 
are consistent with the procedure (A) by the following 
identification 

Ui = ul'^^, for i = 1, 4 

Uj ~ ^^'^ ^ ^ ^' ^' ^' (^^) 

More precisely, we pick the correction terms proportional 
to ^3 in Eq. Ht)2|l . which is estimated as logL in d = 6 — e. 
Putting L = e"*' and takin g th e first order of SI, we get 
the beta functions in Refs.[ll|23. This is an expected 
result, because we can obtain the replicated Hamiltonian 
5'.(,p from the cumulant expansion 




with use of F^qs.^ and ||2Hll- We find that S'^^^ with 
Eq. (gni) gives Eq. 



The procedure (B) contains some diagrams of 0(n). 
These diagrams never appear by way of the present 
method, because those contain a connected component 
without external lines. 



IV. EXPANSION PARAMETERS 
A. Preliminaries 

Since all UjS have the dimension 4 — d, they are ap- 
parently irrelevant when c? > 4. However, coefficients 
of the perturbation series in {uj} are expressed as poly- 
nomial of A having the dimension 2, which can change 
the relevance of ujS. That is, since A is relevant for all 
dimensions, a term containing it can become larger as re- 
peating the RGT. Hence, it is important to explore how 
A appears in the perturbation series. 

To make our argument clear, let us consider St as an 
example. The coefficient of the first order in mi, which 



appears in the first term of St in Eq. 1)62(1 , is depicted by 
the diagrams in FIG. [7| The internal line of (a) in the 
figure carries G{q), which behaves as at the critical- 
ity, while the internal line of (b) brings AG{q)^ ^ A/g^. 
Namely, the dominant contribution is proportional to 




(a) (b) 

FIG. 7: Leading correction from ui to St. 

Aui, as pointed out in the previous literaturesj^ l2lj |. 
We then introduce 

gi=Aui. (67) 

as one of expansion parameters. It has the scaling di- 
mension 6 — d, which implies that the upper critical di- 
mension is six. Therefore we cannot ignore ui in o? < 6, 
even though ui itself is irrelevant when 4 < d. 
In addition, we define 

FICEfa) is then proportional to ui = go<7i- Generally, if 
a diagram includes the irrelevant parameter go, its con- 
tribution is less relevant. 

Next we discuss how A couples with the other param- 
eters U2, U5. Let us look at the first-order terms in U2 
and U3 respectively, which appear in the second term of 
St in Eq. (Note: we have defined U3 = U2 + U3.) 

It turns out that their coefficients do not carry A be- 




U2 us 
(a) (b) 




U2 US 



(c) (d) 

FIG. 8: Diagrams in (a) and (b) give leading correction to t 
by U2 and 113 respectively. However, diagrams in (c) and (d) 
do not contribute to St, because they both have two connected 
components before taking the random-field average. 

cause putting x on the internal line produces discon- 
nected components as shown in FIG. |HIc) and (d), so 
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that (c) and (d) do not appear in the perturbation se- 
ries. Thus no A associates with U2 and u^, so that we 
choose 

33 = A%3 = U3 (69) 

as expansion parameters. They have the scahng dimen- 
sion A — d. 

Let us turn to the parameter u^. It contributes to St 
by the combination uiu^ in the lowest order. See FIG. 
15] Here x is forbidden, for the same reason it is in FIG. 

IHIc). 




FIG. 9: One of the leading contributions from m to St. 

Therefore, this contribution is proportional to uiu^ = 
gi (A~^U4). Then we add 

g4 = A-^U4 (70) 

to the expansion parameters. It has the scaling dimen- 
sion 2 ~ d. 

As for U5, the leading correction to 6t comes from the 
order in ufu^ — gf (A^^us), as we depicted in FIG.ITUl 




FIG. 10: One of the leading contributions from M5 to St. 
Hence we define 

35 = A-2w5, (71) 

which has the dimension —d. 

In this way, we can express the perturbation series of 
dt in terms of go, ...,55. Further, we can rewrite the re- 
cursion equations for A and Ui (1 = 1, 5) into those for 
(7^ (/I = 0, 5). It is easily confirmed from the explicit 
form of Eas. (|59|l and Ht)2|) that no A appears in these 
series at least in the second order in {g^}. Therefore a 
naive dimensional analysis is expected to work. Since the 
scaling dimensions of the new parameters are 

M = -2, [.gi]=6-rf, [52] = [.93] - 4 - rf, 
M = 2-d, [55] = -d (72) 

respectively, only gi is relevant near d = 6. This suggests 
that the extra parameters 32, •■■,55 do not play any im- 
portant role for critical phenomena in the RFIM when d 
is close to 6. 



B. The Gaussian case 

In order to make our argument more precise, we need 
to show that recursion relations for {t, {g^}) do not con- 
tain positive powers of A for all orders. 

To this end, it is instructive to consider the case where 
the probability distribution of impurities is Gaussian, i.e., 
M2, M5 are ignored ,3^.21.]. Suppose that u[ is expressed 
as 

00 

ai — l 

Here fa^ (A, t) is obtained from the sum of all diagrams 
for [w4]p[p] with fli 0^ vertices. Since the number of x in 
a diagram gives the power of A, fa^ (A, t) is a polynomial 
of A with a finite degree. Letting jai be that degree, we 
can write 

/,,(A,t) = ^c„A", (74) 

where c„ is a coefficient that depends on t and ai. Hence 
u'l is written as 

ai—1 n—0 

= L^'^-^'^ J2 ^''"^ < • (75) 

ai—l n—0 

Let US obtain 7^^ . Since the diagrams for u'^ have a 
single connected component, internal lines without x are 
needed at least oi — 1. Excluding four external lines from 
the total lines 4ai, internal lines on which we can put x 
are at most 

i(4ai-4-2(ai-l))=ai-l. (76) 
This is nothing but 71. Thus 

(00 ai — 1 \ 

E5rE^«5o(--^^-" ■ (77) 
ai=l n=0 / 

Note that the power of go, ai — 1 — n, is non- negative in 
the above summation. Similar observations of A' and t' 
lead to the following form: 

Coo ai+l \ 

ai = l n=0 ) 

and 

Coo a\ \ 

^+ Effi'E^^ffo"'"" (79) 
ai=l n=0 / 

where {(i„} and {e„} are i-dependent coefficients. From 
Eqs.ljTTl and GHl, one finds that g'^ = A'tt^, = (A')"^ 
and t' are expanded by g\ and go for all orders. 
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C. General case 



Applying the result of Eq. to Eq. , we get 



Next, we extend the above argument to the case where 
the extra parameters 1*2 j ••■7 ""5 are included. Let Z"*" be 
the set of non-negative integers. Define 

/= {a : a = (ai, as) , flj e Z+, 
5 

Vj = l,...,5, 5]a, >l}. (80) 
i=i 

For a G /, we use the notation 



I!"?- (81) 

The recursion equation for Uj can be written as 



= L4e-3d ^ ^-fi^a g dA^'-^i. (82) 

ael n=0 

Here we compute 7^, the degree of A for m° in u'^. The 
number of connected components in the vertex corre- 
sponding to Uj is denoted by aj, namely, 

ai = 1, a2 = as — 2, = 3, as = 4. (83) 

Here we consider a diagram proportional to u° in Eq. 
(|82|l . Before connecting the vertices in the diagram by 
internal lines, there are 



u'^=L'0-3,^,,Y^gaJ2^^g.i-n^ (89) 

Similarly, we find that A' and t' are expressed by the 
following expansion: 

A' = L'0~'A(l + Y^g^f2dn9r-A 

\ ael n=0 / 



ael n=0 



with 



7a 



(90) 



(91) 



E 



onai 



(84) 



Using Eq. H89|l and the above expansion for A', we con- 
clude that g'j = {A')~^^u'j is expanded by {5^}. Further, 
it is obvious from Eq. (|90|l that and t' also have per- 
turbation series in terms of {5^}. 

It should be noted that some physical quantities, such 
as free energy, are proportional to the relevant parameter 
A. Hence if we compute an exponent associated with 
those quantities, we have to know the singular behavior 
of A near the criticality. As we can see in the first line 
of Eq. (|90|l . A is renormalized multiplicatively and its 
correction is expanded in terms of {.g^}. Therefore, we 
can compute the singular behavior within the framework 
of perturbation in terms of {5^}. 



connected components. Since the diagram has aj con- 
nected components after connecting the vertices, we need 
at least 



E 



i=l 



(85) 



no-x internal lines. Similar computation for Eq. H7t)|) 
leads to 

= ^ ^4^a, -4-2 ^^a.a, -aj^^ 

5 

= - ^ Aa^, 



(86) 



where we have defined (3j = aj — 2. Explicitly, 

/3i = -l, /32-/?3 = 0, /34 = 1, /35 = 2. (87) 

According to the definition of the expansion parameters, 
we find 

gj = A'f^Uj, j = l,...,5. 



D. Observation for dimensional reduction of the 
RFIM 



We have shown that the transformed parameters {g'^} 
can be expressed as positive series in {5^}. Since the 
expansion parameters other than gi are irrelevant near 
d = 6, we may ignore them. 

As we have mentioned in Sec. IIII Dl our recursion re- 
lations are consistent with the replica method studied by 
Brezin and de Dominicis if we adop t the limiting proce- 
dure (A) explained in Sec. IIII Dl 15] . However, it is con- 
cluded that the non-trivial fixed point becomes unstable 
due to M2,...,U5 in Ref . [Tsll . The discrepancy between 
this and our conclusion is resolved as follows: in Ref. [l5| . 
the recursion relations of parameters gj = Auj are com- 
puted. Since [gj] = 6 — d for all j, it is possible that 
some of gj becomes larger as the RGT is repeated. Nev- 
ertheless, it does not mean that a diagram proportional 
to gj brings infrared divergence near the upper critical 
dimension. 

In fact, if we write 6t in Eq. lt)2|) by gj instead of Uj, 
one can easily check that gjS for j = 2, 5 are always 
combined with A~^ or A~^. Thus, even though gj (j = 
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2, ...,5) behave as g'^ ^ L'^^^gj with Cj > 0, the negative 
powers of A suppresses growth of terms proportional to 
(jj in St. The discussion in the previous subsection shows 
that the association with the negative powers of A occurs 
for all orders. Thus gjs {j = 2,..., 5) do not contribute 
to exponents. 

Although the extra parameters remain irrelevant near 
the upper critical dimension, it is plausible that those 
parameters become relevant when e exceeds some finite 
value Cc, which can cause the breakdown of the dimen- 
sional reduction in c? = 3. The existence of such critical 
value is consistent with a high-temperature expansion by 
Houghton et al^^. They concluded that dimensional 
reduction occurs in c? = 5 and 6 while the phase tran- 
sition becomes first order in d < 4. It is strongly sug- 
gested that 1 < ec < 2. On the other hand, another 
high-temperature expansion performed by Gofman et al 
suggests that the breakdown of dimensional reduction in 
d < 5(12i]- It can be interpreted as < ec < 1. 

As we explained above, dimensional reduction can sur- 
vive for sufficiently small e. In this case, the exponents 
J/, ry and fj are calculated as|25| 

J^=i + ^e + 0(e2), rj^fj^j^e^ + Oie^) (92) 

However, other parameters not considered here may 
bring the breakdown of dimensional reduction for all 
d < 6. It is conjectured by Feldman that the ap- 
parently higher-dimensional operators '^^bi^'^ ~ 0'')^' 
{I > 1/e^) turn to relevant onesp^. In the Harris- 
Lubensky method, similar operators X]ab('^°'^'')' i'^" 
troduced, taking into account the following random av- 
erage: 

[vi (fci, ...,ki);vi {k[, ...,/c,)]p[^] = tl;,5 (fci H h k'l) . 

(93) 

The parameter ui is, itself, irrelevant. However, follow- 
ing the argument in Sec. IIVI one can find an expansion 
parameter proportional to ui is (A)' ^u;, which has the 
canonical dimension {I — 1)(4 — d) ^ 1{A — d) for large I. 
Feldman shows that they acquire anomalous dimension 
with O (?^e^) in the second-order perturbation. It means 

that (A)' ^ ui can transmute a relevant parameter for 
sufficiently large I satisfiying I > 1/e^. The conjecture 
should be checked by a non-perturbative means that is 
beyond the scope of this paper. Nevertheless it is consis- 
tent with our result in that the operator with I = 2 (i.e., 

^ ((^"0'') ) is irrelevant for sufficiently small e. 

V. SUMMARY AND DISCUSSION 

In this paper, we have studied the critical phenomena 
of the extended Ginzburg-Landau theory for the random- 
field Ising model. We have employed the renormaliza- 
tion group for the probability distribution of the impu- 
rities. Probability distribution is characterized by non- 
trivial cumulants that bring extra parameters, which are 



essentially identical to the new coupling constants in the 
replica Hamiltonian introduced in Refs.fTslIT^. In con- 
trast to the replica method, our approach does not re- 
quire any limiting procedures, and hence no artificial am- 
biguities arise. Thus we can definitely determine the scal- 
ing dimensions of expansion parameters. We have found 
that extra expansion parameters do not affect the critical 
phenomena in c? = 6 — e with sufficiently small e. On the 
other hand, those parameters could be an obstruction to 
the dimensional reduction at some finite e. This result 
indicates that we cannot rule out dimensional reduction 
near d = 6 by including the extra coupling constants 
in Eq. (|28|l . It is consistent with the high-temperature 
expansion by Haughton et aZ[3|. On the other hand, 
another high-temperature expansion by Gofman et al is 
not consistent with Ref. 16], which may indicate that the 
dimensional reduction does not occur in any d < 6|l2j|. 
It is important to resolve this discrepancy for clarifying 
mechanisms in the phase transition in the RFIM near the 
upper critical dimension. 
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APPENDIX A: COMPUTATION OF THE 
RECURSION EQUATIONS 

In this appendix, we show technical details of comput- 
ing the correction terms Eq. (|59|l in O {uiUj). Eq. H()2|) 
can be obtained in a similar manner. 

We do not obtain an explicit form of Vj. Instead, we 
derive all possible diagrams in the desired order that ap- 
pear in the right-hand side of FIG. |S1 

First, we consider diagrams without x on an internal 
line. As we have presented in FIG. 13 the cumulants can 
be expanded by the five kinds of vertices. The lowest- 
order correction to u'^ comes from the vertices with two 
external legs. There are nine such vertices as depicted in 
FIG. ^2 Merging two of them, we obtain diagrams that 
contribute to M^(j = 1,...,5) up to the order 0{uiUj). 
We write the number j in Table if the resulting dia- 
gram contributes to u'^. For instance, we can read from 
Table m that the diagram made of the vertices A and B 
contributes to 1*2- As we have mentioned in Sec. IIII Al 
any connected components of a diagram in the right-hand 
side of FIG.Elmust contain at least one external line. In 
Table |l| a blank means that the corresponding diagram 
does not satisfy this condition. 

Now we calculate the numerical factor to each diagram, 
as we have demonstrated in llll Bl In contrast to the pure 
(/)^ theory, a one-loop diagram has multiple connected 
components in general. There is the combinatoric factor 
associated with each connected component, which can be 
computed in the same way as with the pure theory. Here 
we denote the product of the combinatoric factor of each 
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X X 



B 



C 



D 



FIG. 11: nine vertices for u'j {j = 1, ...,5) 



ABCDEFGHI 



A 1 

B * 

C * 

D * 

E * 

p * 

G * 

H * 

/ * 



3 4 

4 5 
4 5 



TABLE I: Admissible one-loop diagrams made from the ver- 
tices in FIG, mi If a diagram contributes to u'j, the number 
j is entered. 



connected component as 



(Al) 



Next, consider the case where a cumulant for Uj has 
two (or more) identical potential Vi. Suppose that v[ 

is expressed by the series v^^^^ . If a diagram for m' 



consists of and v)^ ' with k ^ k', then 



(fc') 



\p[p] 



-'i ' ■ 



„(fc'). 



„(fc'). 



,(fc). 



Such multiplicity is denoted as 

nR. 



Plp] 
(A2) 

(A3) 



Finally, we have to take into account the number of ways 
of contraction in the random-potential average. Here it 
is denoted as 

nc- (A4) 

Then, the factor with a diagram n^? is computed as 

np = UAriRnc- (A5) 

In Tables llil VII we outline those factors. The first row 
denotes the power of A of the diagrams. The second 
row shows the name of the vertices in FIG. that are 
ingredients of the diagrams. The last row indicates the 
parameter dependence of the diagrams. 





A" 


A' 




AA AC AD 


AA 


riA 
riR 

nc 
np 


-i 6 6 
111 
111 

-I 6 6 


-3 
1 
1 

-3 




-2 


UiU2 



TABLE II: Diagrams for u'l and their numerical factors 





A" 


A' 




AB 


AE 


BC 


BD 


CD 


AC AB 




3 
2 


3 


3 


3 


3 


-3 -3 


UR 


1 


1 


1 


1 


1 


1 1 


nc 


1 


1 


1 


1 


1 


1 1 


np 


3 
2 


3 


3 


3 


3 


-3 -3 




U1U2 


U1U4, 


U2U3 


u'i 


U2U3 


U1U3 U1U2 



TABLE III: Diagrams for u'2 and their numerical factors 





A" 


A 


i 


A^ 




AF 


AH 


CC DD 


CF 


DF 


AD 


AF 


AA 




1 

2 


i 
2 


1 


1 


1 


1 


-1 


-1 


1 
4 


Ur 


2 


2 


1 


1 


2 


2 


1 


2 


1 


nc 


1 


1 


2 


1 


1 


1 


2 


1 


2 


np 


-1 


1 


2 


1 


2 


2 


-2 


-2 






U1U3 


U1U4 


ui 


ui 


ui 


U2M3 


U1U2 


U1U3 





TABLE IV: Diagrams for U3 and their numerical factors 

It is worthwhile mentioning the relationship to the 
replica method J^]. From Sl^^ defined in Eq. (|64|l . 



dcj)ad(j)i3 



rep 



rep 



'■2 
rep 



2 

rep 



rep 



fiiy 
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A" 




AG 


AI 


BB 




BF 


BH 


C£ 


CG 


Ci/ 


riA 


i 

2 


i 

2 


1 

2 


2 


i 
2 


i 
2 


1 


1 


1 


riR 


1 


1 


1 


1 


2 


2 


2 


1 


1 


nc 


1 


1 


2 


2 


1 


1 


2 


1 


1 


np 


1 
2 


1 
2 


-1 


4 


-1 


1 


4 


1 


1 






M1M5 


M2 


U2«4 


U2U3 


U2U4 


U3M4 


M3M4 


U3M4 





A" 






A 


i 






DE 


DG 


EF 


AE 


BG 


BD 


GD 


GG GF 




1 


1 


1 


-1 


-2 


-1 


-1 


-1 -1 


riR 


2 


1 


2 


2 


2 


2 


2 


1 2 


nc 


1 


1 


1 


1 


1 


1 


1 


2 1 


np 


2 


1 


2 


-2 


-4 


-2 


-2 


-2 -2 




M2W4 


1121*4 


U3«4 


U1U4, 


U2U3 


«^ 


M2'"3 


2 2 
M3 1^3 





A^ 


A^ 




AG 


BB 


BF 


AC AB 




~1 


-1 


-1 


1 1 

2 4 


nii 


1 


1 


2 


1 2 


nc 


1 


2 


1 


2 2 




-1 


-2 


-2 


1 1 




Ullt4 


"2 


M2W3 


111^3 U1U2 



TABLE V: Diagrams for 1*4 and their numerical factors 





A" 


A^ 




BG 


B7 


GI 


F^ 


SG 


BE 


GG GF 


BG 


riA 


1 

2 


i 
2 


1 


1 


1 


-1 


-2 -i 
^ 2 


-1 


riR 


4 


4 


6 


6 


4 


12 


12 12 


4 


nc 


3 


3 


1 


2 


3 


2 


1 2 


3 


np 


-6 


6 


6 


12 


12 


-24 


-24 -12 


-12 




W2lt4 


U2U5 


U3U5 


U4 




112M4 


ti3lt4 U3114 


M2W4 





A^ 




BG 


GG BB 


UA 


i 
2 


1 k 


nR 


12 


1 6 


nc 


2 


6 4 


np 


12 


6 6 




M2M3 


2 2 

M3 U2 



TABLE VI: Diagrams for u'^ and their numerical factors 



These terms exactly produce the vertices in FIG. ^2 un- 
der the condition that each line be connected in the case 
where the corresponding field has the common repUca 
index. The one-loop diagrams identical with Tableware 
obtained using the replica propagator 

We can check that the recursion equations identical with 
those in the main text are derived by this method, if we 
take the limit n ^ with UjS fixed |2^. 



APPENDIX B: CRITICAL EXPONENTS NEAR 
THE UPPER CRITICAL DIMENSION 

In this appendix, we show an outline of computing 
critical exponents in d = 6 — e. 

Let us introduce the connected and the disconnected 
two-point function gdk] /i) and gd{k] /i) by the following 
formulae: 

[((/>(fci) ;(/)(fc2))]p[p] = gc {ki; fi) S {ki + ^2) 

[(</) (fci)} ; (</) (fc2))]p[p] = gd (fci; m) S (fci + fc2) ,(B1) 

where fi represents a point in the parameter space /i — 
{t, go, g^). Since gdk;^) is computed from the sum 
of connected diagrams that appear in St, it can be ex- 
panded by {5^}. On the other hand, gd (fc; /i) is made of 
diagrams having two connected components appearing in 
SA. Namely, gdik; ^) can be written as 

gd (fc; = A.9d (fc; A^) , (B2) 

where gd (k\ /i) is expressed as a perturbative series of 
{5^}. Suppose that A transforms as 

A' 2± L'^A. (B3) 

The correlation functions satisfy the following transfor- 
mation law: 

gap;M) = i''''.9c(^p;M') 

A~gd {p; m) = L^'-'+'^A-gd {Lp; m') • (B4) 

First, we focus on the critical exponents rj and fj, which 
determine the small-momentum behavior of correlation 
functions at the criticality: 

9c{p\p) ^ {pi +P2) 

9d{k;fi) ~ -^Hpi+P2)- (B5) 
p^ 

Suppose that we apply the RGT n times at the criti- 
cality, where n satisfies 

L>i=Aei, (B6) 

with ei being some d dimensional unit vector. The prob- 
ability distribution P approaches the fixed-point distri- 
bution, P*, characterized by /i*. We can evaluate gdPy 1^) 
as 

9c{pi\p) 

^ (^-j ffc(Aei;A**), (B7) 

where /i*-"-* specifies the probability distribution of impu- 
rities after having applied the RGT n times. Comparing 
the definition (jB5|) . we have 

e^\{2 + d~i^). (B8) 
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Applying the same argument to gdip] m)j obtain 



and 



= -{4 + d-K-f]) 



K — 2 + 7] — rj. 



(B9) 



(BIO) 



The exponents 7 and 7 are computed, respectively, 
from gc{0; fx) and 5^(0; Ij) in the disordered phase. In this 
case, the RGT is repeated n times, where L" is equal to 
the correlation length ^. Then 



X = 5c(0;Ai) 



(Bll) 
(B12) 
(B13) 
(B14) 



Assuming that 

we get, with the help of Eq. (|B8p . 
where 

7 = (2 - 77)1/. 
Similar arguments for (7d(0; /i) lead to 

X = <?d(0;/i)^(t-tc)-^ (B15) 

with 

^^{A-ff)v. (B16) 

Next we consider the exponent a. The singular part 
of free energy transforms as 

F{y) = F{pi'), (B17) 

hence its density / transforms as 

/(mO - ^ - = ^'^/(A^)^ (B18) 

where /(yit) has the form of 

f{^l) = A/(m) (B19) 
with /(/i) having the perturbative series of {5^1}. Thus 

f{^l) ~ - r'^^", (B20) 

which means that 

2 - a = (d - 2 - r/ + ^) . (B21) 

Eq. l|B20p shows that n is identical with the exponent 
of the singular part of the free energy in a correlation 
volume which is often denoted by 

In d = 6 — e, we can perform the perturbation 
explicitlydHll^ . We begin with the equality 

{kl + t') 5 {ki + k2) = [V2(fcl,fc2)]p,[p] 

= K(fci,fc2)]p[p]. (B22) 



Employing Eq. (|58|l and denoting 

[(5w2 (pi,P2)]p[p] =(5r2(pi;/i)(5(pi+p2), (B23) 
we get 

W2 (fci,fe)]p[p] = 

L^e-'i (fc2^-2 +t + ST2 (pi; A*)) 5 (fci + fca) (B24) 



Comparing the coefficient of k\ in Eas. (|B22|l and HB24|I . 
we get 



= 1 



d 
dpi 



5^2 (pl;fJ-* 



(B25) 



Pi=0 



which determines the exponent 77. On the other hand, 77 
is computed by a correction to A. Define STi^{p; fi) by 

K {h);v[ (fc2)]p[p] = 

^2ie-d)^ (1 + ^Fa (pi; m)) -5 (pi +P2) • (B26) 



We can readily derive 

A' ^L^^^'^ A {I + ST A (0;/i*)) 



(B27) 



Repeating similar calculations for deriving Eq. HB25() , we 
obtain 



= 1 + (SEa (0;/i*). 



(B28) 



When d = 6 — e, we can evaluate dT2 and ST a by the e 
expansion. Here we can ignore the irrelevant parameters 
(/2, ...jffs. In this case, we have 



d 
dpi 



Pl=0 



ST2{pl;fi)=STA (0;/i) 



(B29) 



where // is a point in the parameter space satisfying 172 = 
• • • = (75 = 0. The above equation indicates that 

= 77 (B30) 

with Eqs. (|B25|) and HB28I) . This equality shows that 

K = 2 (B31) 

from Eq. (jBlOl) . It derives 



9'1-L 



K+4-d-2?7 



(51 +5gi) 



-6-(i-277 



(.91 + 5gi). 



(B32) 

where Sgi is shown to be identical with that of the 4 — e 
dimensional pure 0* theory, with the coupling constant 
Further, the correction term St is also equal to that 
of the pure theory in 4 — e dimensions. In this way, we 
can rederive the result of dimensional reduction near the 
upper critical dimension. 
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